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ABSTRACT 

We present new rates for the 22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg reactions, with 
uncertainties that have been considerably reduced compared to previous estimates, and 
we study how these new rates affect the production of the heavy magnesium isotopes in 
models of intermediate mass Asymptotic Giant Branch (AGB) stars of different initial 
compositions. All the models have deep third dredge-up, hot bottom burning and mass 
loss. Calculations have been performed using the two most commonly used estimates 
of the 22 Ne + a rates as well as the new recommended rates, and with combinations of 
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their upper and lower limits. The main result of the present study is that with the new 
rates, uncertainties on the production of isotopes from Mg to P coming from the 22 Ne 
+a-capture rates have been considerably reduced. We have therefore removed one of 
the important sources of uncertainty to effect models of AGB stars. We have studied the 
effects of varying the mass-loss rate on nucleosynthesis and discuss other uncertainties 
related to the physics employed in the computation of stellar structure, such as the 
modeling of convection, the inclusion of a partial mixing zone and the definition of 
convective borders. These uncertainties are found to be much larger than those coming 
from 22 Ne + a-capture rates, when using our new estimates. Much effort is needed to 
improve the situation for AGB models. 

Subject headings: stars: AGB and post-AGB — nuclear reactions, nucleosynthesis, 
abundances 

1. Introduction 

The origin of the stable magnesium isotopes, 24 Mg, 25 Mg and 26 Mg, is of particular interest to 
astrophysics because Mg is one of the few elements for which we can obtain isotopic information 
from stellar spectroscopy. The ratio 24 Mg: 25 Mg: 26 Mg has been derived from high-resolution spectra 
of cool dwarfs and giants in the thin and thick disk of the Galaxy (Gay & Lambert 2000; Yong, 
Lambert & Ivans 2003b), and for giants stars in the globular cluster (GC) NGC 6752 (Yong et 
al. 2003a). These observations show that many of the stars, including relatively metal-poor stars 
([Fe/H] < —1-0), have non-solar Mg isotopic ratios 1 with enhancements in the neutron-rich isotopes, 
25 Mg and 26 Mg, compared to what is expected from galactic chemical evolution (GCE) models. The 
main stellar nucleosynthesis site for all three stable isotopes is hydrostatic burning in the carbon and 
neon shells of massive stars that explode as Type II supernovae (Woosley & Weaver 1995; Chieffi &: 
Limongi 2004) . The abundances of the neutron-rich Mg isotopes are further enhanced by secondary 
a-capture processes operating in the helium shell. The amount of 24 Mg produced does not strongly 
depend on the initial metallicity of the model and is an example of primary nucleosynthesis 2 , 
whereas the amounts of 25 Mg and 26 Mg produced scale with the initial metallicity of the star. This 
means that the Mg content of the ejecta of low metallicity supernovae will be mostly 24 Mg with 
very little 25 Mg and 26 Mg produced, with typical ratios 24 Mg: 25 Mg: 26 Mg « 99.0:0.50:0.50 from a 
25M Q supernova model with metallicity Z = 0.01Z Q (Woosley & Weaver 1995). 

Previous studies have shown that GCE models using ejecta from massive stars match the obser- 



ve solar Mg isotopic ratios are 24 Mg: 25 Mg: 26 Mg = 79:10:11 (Lodders 2003). 

2 Primary production means that the species is produced from the hydrogen and helium initially present in the 
star and the amount produced is relatively independent of the metallicity, Z, whereas secondary production requires 
some heavier seed nuclei to be present, and the amount produced scales with Z. 
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vational data well for [Fe/H] > — 1.0 but severely underestimate 25 ' 26 Mg/ 24 Mg at lower metallicities 
(Timmes, Woosley &; Weaver 1995), indicating that another production site for the neutron-rich 
Mg isotopes at low metallicities is required to account for the observations. Recently, Fenner et 
al. (2003) included the predicted stellar yields of 25 Mg and 26 Mg from Asymptotic Giant Branch 
(AGB) stars (Karakas & Lattanzio 2003) along with yields from Type II supernovae (Woosley &: 
Weaver 1995; Limongi & Chieffi 2002) into a GCE model of the solar neighborhood. The GCE 
model with the AGB contribution could successfully match the Mg isotopic ratios of the metal-poor 
Galactic disk stars while the model without an AGB contribution could not. This result indicates 
that low-metallicity intermediate mass AGB stars may play an important role in the production of 
these species in galaxies and stellar systems. The production of the neutron-rich isotopes in AGB 
stars is also of interest in relation to the non-solar Mg isotopic ratios observed in giant stars in 
globular clusters (Yong et al. 2003a; Yong, Aoki & Lambert 2005). The non-solar Mg isotopic ratios 
observed in NGC 6752 have been attributed to AGB stars, but Fenner et al. (2004) used a GCE 
model with tailor made AGB yields from Campbell et al. (2005) and failed to match the abundance 
patterns observed in stars in that cluster. As Ventura & D'Antona (2005a,b) have pointed out, 
there are still many uncertainties that effect the stellar yields, so an AGB solution to the globular 
cluster anomalies cannot be ruled out at present. 

Further motivation for the study of the production of the Mg isotopes in AGB stars is given 
by their relevance in the important current debate on the apparent variation of the fine-structure 
constant (Murphy et al. 2001; Ashenfelter, Mathews & Olive 2004; Fenner, Murphy & Gibson 
2005), and the origin of pre-solar spinel grains, some of which show enhancements in both 25 Mg 
and 26 Mg compared to solar (Zinner et al. 2005). 

Briefly, intermediate- mass stars (initial mass ~4 to 8M ) will enter the thermally-pulsing phase 
with a hydrogen (H)-exhausted core mass (hereafter core mass) M c > 0.8M Q , after experiencing 
the second dredge-up (SDU) (Lattanzio et al. 1996; Busso et al. 1999; Herwig 2005). The SDU 
brings the products of H-burning to the stellar surface (mostly 4 He and 14 N), and will slightly 
alter the composition of the Mg isotopic ratios, with an enrichment in 26 Mg at the expense of 
25 Mg. For the massive, Z = 0.0001 models, the SDU is the first time the surface abundances 
are altered, because there is essentially no first giant branch phase (Herwig 2004a). Following the 
SDU, the He-burning shell becomes thermally unstable and flashes every few thousand years or 
so. The energy from the thermal pulse (TP) drives a convective pocket in the He-rich intershell, 
which thoroughly mixes the products of He-nucleosynthesis within this region. Following the TP, 
the convective envelope moves inward in mass and may reach the region previously mixed by the 
flash-driven convective pocket. This mixing event is known as the third dredge- up (TDU) and if it 
occurs, is responsible for enriching the envelope in material from the H-exhausted core. Following 
the TDU the star contracts and the H-shell is re-ignited and provides nearly all of the surface 
luminosity for the next interpulse period. The thermal pulse - TDU - interpulse cycle may occur 
many times during the TP- AGB phase; how many times depends on a number of factors including 
the convective model which determines the surface luminosity and mass-loss rate and hence the 
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total AGB lifetime (Ventura & D'Antona 2005a). 

Hot bottom burning can also occur when the base of the convective envelope becomes hot 
enough to sustain proton-capture nucleosynthesis (Lattanzio et al. 1996). If temperature at the 
base of the envelope is sufficiently hot (over ~ 60 x 10 6 K), the NeNa and MgAl chains may operate 
alongside the CNO cycle; Li production is also possible via the Cameron- Fowler mechanism (Sack- 
mann & Boothroyd 1992; Lattanzio et al. 1996), which operates at lower temperatures (typically 
~ 30 — 40 x 10 6 K). Frost et al. (1998) noted that intermediate-mass AGB stars may become lumi- 
nous, optically obscured carbon stars near the end of the TP-AGB, when mass loss has removed 
much of the envelope, extinguishing HBB but allowing dredge-up to continue. 

Karakas & Lattanzio (2003) described in detail the various nucleosynthesis processes that alter 
the Mg isotopic ratios in AGB stars. To summarize, 25 Mg and 26 Mg are synthesized in the He-shell 
during thermal pulses by the reactions 22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg, when the temperature 
exceeds about 300 xlO 6 K. The amount of Mg produced depends on the thermodynamic conditions 
inside the pulse as well as on the composition of the intershell, which will have been altered by 
previous H and He-burning. Neutron captures, in particular the 25 Mg(n, 7) 26 Mg reaction, can also 
alter the Mg isotopic ratio in the intershell, where the neutrons come from the 22 Ne(a, n) 25 Mg 
reaction (Herwig 2004b). HBB can also significantly alter the surface Mg isotopic ratio via the 
activation of the MgAl chain, which can result in the destruction of 24 Mg if the temperature 
exceeds ~ 90 x 10 6 K. 

The stellar yields of 25 ' 26 Mg presented in Karakas Sz Lattanzio (2003) and shown in figure 1 
were calculated from models covering a range in mass (1M to 6M Q ) and metallicity ([Fe/H] 
= 0,-0.3,-0.7). From this figure we see that the most massive AGB models produce the most 
25 > 26 Mg, as a consequence of higher temperatures in the He-shells compared to lower mass stars; 
we also notice an increase in production at a given mass with a decrease in metallicity. The 
computationally demanding nature of AGB models precluded a detailed study in that paper of 
the effect of the major uncertainties (mass loss, nuclear reaction rates, convection). A recent 
comprehensive study by Ventura Sz D'Antona (2005a, b) demonstrated that the predictive power 
of AGB models is still seriously undermined by these uncertainties. The theory of convection 
has a significant effect on the structure and nucleosynthesis (Ventura Sz D'Antona 2005a), whilst 
varying the mass-loss rate results in larger changes to the stellar yields than varying the nuclear 
reaction rates (Ventura Sz D'Antona 2005b). However, the magnitude of the errors associated 
with the relevant nuclear reaction rates is still one of the key questions concerning predictions of 
magnesium production in intermediate-mass AGB models. Whilst the rates of the 14 N(a,7) 18 F 
and 18 0(a,7) 22 Ne reactions are well determined (Gorres et al. 2000; Dababneh et al. 2003), the two 
key a-capture reaction rates: 22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg, suffer from large uncertainties 
at the stellar energies appropriate for AGB stars (Koehler 2002). For example, at typical He-shell 
burning temperatures, T « 300 x 10 6 K, the NACRE compilation (Angulo et al. 1999) give an upper 
limit to the 22 Ne(a, n) 25 Mg reaction rate that is about 47 times larger than the recommended rate. 
At lower temperatures, the uncertainties are even larger. 
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The aims of this paper are two fold. First, we present new reaction rates for the two key a- 
capture reactions, with considerably reduced uncertainties compared to those given by the NACRE 
compilation. Second, we use these new rates within models of different metallicities ([Fe/H] ~ 
0, —0.3, —0.7, —2.3), for a typical mass (5M ) that produces the Mg isotopes during the thermally- 
pulsing AGB (TP-AGB) phase. For each model we calculate the stellar yields and compare to 
previous nucleosynthesis calculations using older estimates of these reaction rates, including those 
by Herwig (2004b) and Ventura & D'Antona (2005a,b). We also examine the effect of other model 
uncertainties, in particular the inclusion of a partial mixing zone and mass loss, on the stellar yields. 

The paper is organized as follows. In §2 we discuss the numerical method used for the stellar 
model calculations, including a discussion of the input physics used. In §3 we present new rates for 
the 22 Ne +a capture reactions. The results of the calculations are presented in §4 and discussed in 
55. 



2. The numerical method 

We calculate the structure first and perform detailed nucleosynthesis calculations afterward, 
using a post-processing algorithm. The stellar structure models were calculated with the Monash 
version of the Mount Stromlo Stellar Structure Program; see Frost & Lattanzio (1996) and refer- 
ences therein for details. Mass loss on the first giant branch is included using the Reimer's mass-loss 
prescription (Reimers 1975) with the parameter r] = 0.4; on the AGB we use the formulation given 
by Vassiliadis and Wood (1993) in all models unless indicated otherwise. We calculate two models 
(5M Q , Z = 0.02 and Z = 0.0001) using Reimer's mass loss on the AGB with the parameter r] = 3.5. 
Here we assume that the Vassiliadis and Wood (1993) and Reimers (1975) mass-loss prescriptions, 
derived for solar- like metallicities or those appropriate for the Magellanic Clouds, can be applied 
to Z = 0.0001 intermediate-mass AGB models without modification. The use of either mass-loss 
law results in a large number of thermal pulses (up to ~ 140). 

All models were calculated from the zero-age main sequence to near the end of the TP-AGB 
phase. The occurrence of the third dredge up and hot bottom burning (HBB) depend critically 
upon the convection model used (Frost & Lattanzio 1996; Mowlavi 1999; Ventura & D'Antona 
2005a) and the method for determining convective borders. Briefly, we use the standard mixing- 
length theory (MLT) for convective regions, with a mixing-length parameter a = I /Hp = 1.75, and 
determine the border by applying the Schwarzschild criterion. Hence we do not include convective 
overshoot, in the usual sense. We do however, search for a neutral border to the convective zone, 
in the manner described in detail by Lattanzio (1986); Frost & Lattanzio (1996) and Karakas et al. 
(2002). We note that this method has been shown to increase the efficiency of the third dredge- up 
compared to models that strictly use the Schwarzschild criterion. 

We performed detailed nucleosynthesis calculations using a post-processing code which includes 
74 species and time-dependent diffusive mixing in all convective zones (Cannon 1993). The details 
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of the nucleosynthesis network are outlined in Lugaro et al. (2004), but we remind the reader 
that we include 59 light nuclei and 14 iron-group species. We also add the fictional particle g, to 
count the number of neutron captures occurring beyond 62 Ni (Lattanzio et al. 1996; Lugaro et al. 
2004). Initial abundances are taken from Anders & Grevesse (1989) for the Z = 0.02 models and 
we assume scaled-solar for the Z = 0.0001 models. We assume an initially a-enhanced mixture 
typical of thin-disk stars (Reddy et al. 2003) for the 5M , Z = 0.008,0.004 models. The initial 
Mg isotopic ratios (e.g. 25 Mg/ 24 Mg) are 21.9% and 35.2% less than solar in the Z = 0.008 and 
Z = 0.004 models, respectively. 

The bulk of the 506 reaction rates are from the REACLIB data tables (Thielemann et al. 1986), 
based on the 1991 updated version. Some of the proton, a and neutron capture reaction rates have 
been updated according to the latest experimental results, see Lugaro et al. (2004) for details, and 
we calculate one set of models using these rates. We define this set to be our standard set and the 
results from these computations are shown in figure 1. In all the other sets of computed models 
we have updated the proton capture rates for the NeNa and MgAl chains to those recommended 
by NACRE (see table 7 in §A for details of which proton capture rates that have been updated). 
With this new set of updated rates, we study the effect of changing the 22 Ne +a rates: we run 
models with our standard choice for these rates (Kappeler et al. 1994, without the inclusion of the 
633 KeV resonance, indicated as K94 in table 6), with the NACRE recommended rates (indicated 
as NACRE in table 6), and with our new recommended rates, as well as using combinations of the 
new upper and lower limits. 

3. Nuclear structure of 26 Mg and the uncertainties in the 22 Ne + a reaction rates. 

The reaction rates for the 22 Ne+a fusion processes are determined by the level structure of the 
compound nucleus 26 Mg above the a threshold at T Q =10.615 MeV and near the neutron threshold 
T n =11.093 MeV. There is only very limited experimental and theoretical information available 
about possible natural parity resonances in that energy range and both the 22 Ne(a, 7) 26 Mg as well 
as the 22 Ne(a, n) 25 Mg reaction rates have substantial uncertainties in the temperature range of 
stellar helium burning. These uncertainties have been a matter of debate for quite some time. 

The uncertainty for the 22 Ne(a, n) 25 Mg reaction is mainly determined by the possible contri- 
bution of a low energy 1~ resonance at 0.538 MeV which has been observed both in photon induced 
neutron emission (Berman et al. 1969) as well as in neutron capture measurements (Weigmann et 
al. 1976; Koehler 2002). Yet it is not clear how strongly this level is populated in the a channel. 
Direct 22 Ne(a, n) 25 Mg measurements (Harms et al. 1991; Drotleff et al. 1993; Jaeger et al. 2001) 
were performed over the entire energy range down to the neutron threshold, but these measure- 
ments were handicapped by cosmic ray induced background in the neutron detectors, therefore 
only an upper limit for the resonance strength was obtained. Also in a 22 Ne( 6 Li,d) 26 Mg a transfer 
measurement (Giesen et al. 1994) only an upper limit for the a spectroscopic factor S a <0.2 was 
determined. The NACRE compilation (Angulo et al. 1999) has based the upper limit for the reac- 
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tion rate on the direct 22 Ne(a, n) measurement by Drotleff et al. (1993). For the lower limit any 
contribution of this resonance was neglected. While more recent 22 Ne(a, n) studies have reduced 
the upper limit substantially (Jaeger et al. 2001) the proposed value is still substantially higher 
than suggested by the a transfer data. A recent systematic re-analysis by Koehler (2002) of the 
25 Mg(n, 7) measurement by Weigmann et al. (1976) pointed out on the basis of a careful R-matrix 
analysis of 25 Mg(n, 7) 26 Mg data that at least three or more additional natural parity states could 
be expected above the neutron threshold. It was shown by Koehler (2002) that in particular the 
possible contribution of the 2 + level at 11.112 MeV to the reaction rate does raise its uncertainty 
by more than an order of magnitude. In view of this debate and in absence of new data it seems 
worthwhile to review the present status and suggest more reliable limits for the rate taking into 
account all of the presently available experimental information. 



3.1. The reaction rate of 22 Ne(a, n) 25 Mg 

As it has been outlined by Koehler (2002), several natural parity resonances near the neutron 
threshold can contribute significantly to the 22 Ne(a, n) 25 Mg reaction rate increasing the neutron 
production at low energies. The upper limits given by Koehler (2002) were scaled to the experimen- 
tal upper limits of the resonance state at £^ m =0.538 MeV (£^=11. 153 MeV, J n =l~) by Jaeger et 
al. (2001). As pointed out above a transfer measurements suggest a substantially lower a strength 
for the £'^=11.153 MeV level. While the experimental lower limit indicates a resonance strength 
of ujj <60 neV the transfer data reduce the upper limit considerably to ujj <27 neV (see also 
Kappeler et al. 1994) after renormalization of the reference resonance. 

Based on the data by Giesen et al. (1994) and Koehler (2002) we have reanalyzed the strength 
of the low energy resonances and derived the resonance strength from the a spectroscopic factors 
normalized to the newly determined strength of the well known resonance at £ , „ m =0.713 MeV 
(-E x =11.328 MeV, J 7V =1~) (Jaeger et al. 2001). The procedure has been outlined in Giesen et al. 
(1994). Table 1 summarizes the level parameters adopted for determining the strengths as well as 
the upper and lower limits for the neutron unbound states between 11.1 and 11.4 MeV in 26 Mg. 
While neutron capture data (Weigmann et al. 1976; Koehler 2002) indicate a large number of natu- 
ral parity states between 11.16 and 11.3 MeV excitation energy, no appreciable a strength has been 
observed for these levels in the 22 Ne( 6 Li,d) 26 Mg transfer reaction. For these states levels we adopt 
an average spectroscopic strength of C 2 S Q =0.005. For the two levels at 11.113 and 11.153 MeV 
at the neutron threshold the experimental data indicate some a strength, which however could 
not be reliably determined due to the lack of experimental resolution in the deuteron spectrum. 
We extracted S a spectroscopic factors of 0.005 and 0.01, respectively with an upper limit of 0.01 
and 0.02. For the neutron and 7-partial widths we adopted the values given by Koehler (2002) 
for calculating the resonance strengths. Using these resonance strengths listed in table 2 we have 
calculated the reaction rate contributions of the resonances for 22 Ne(a, n) 25 Mg. 
The reaction rate is shown as a function of temperature in figure 2. In the temperature range 
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between 0.3 and 0.5 GK the reaction rate is dominated by the contribution of the lowest measured 
resonance at E^ m =0.713 MeV (Jaeger et al. 2001). At higher temperatures the rate is determined 
by higher energy resonance contributions resonances (Harms et al. 1991; Drotleff et al. 1993). The 
threshold resonance states determine the rate at lower temperatures < 0.3 GK. The largest uncer- 
tainty is obviously associated with this low temperature range. Based on the presently available 
data the reaction rate is uncertain by roughly one order of magnitude. This uncertainty depends 
mainly on the uncertainty in the spectroscopic a strengths of the threshold resonances. More 
detailed experimental work is necessary to reduce the present uncertainties. 



3.2. The reaction rate of 22 Ne(a, 7) 26 Mg 

It has been argued before (Kappeler et al. 1994) that the 22 Ne(a, 7) 26 Mg reaction may compete 
strongly at low temperature with the 22 Ne(a, n) 25 Mg capture process. This possibility depends 
sensitively on the contribution of resonances at energies below the neutron threshold. Koehler 
has argued that the strong resonance at £'^ m =0.711 MeV (£^,=11.326 MeV, J 7r =l~), which was 
observed by Wolke et al. (1989) does not correspond to the resonance level observed in the competing 
22 Ne(a,n) reaction channel (Harms et al. 1991; Drotleff et al. 1993; Giesen et al. 1994; Jaeger et al. 
2001). Only resonance levels with small neutron widths T n < T 7 can have an appreciable resonance 
strength in the 22 Ne(a, 7) channel. 

We have recalculated the resonance strengths ujj of all possible resonances in the 22 Ne(a, 7) 
channel. The resonance strengths scales directly with the spectroscopic factors S a (Giesen et al. 
1994), which have been re-normalized to the experimentally determined strength of the state at 
11.328 MeV (Jaeger et al. 2001). The spin parity assignment for the low energy sub-neutron 
threshold resonances is based on the DWBA analysis of the 22 Ne( 6 Li,d) 26 Mg data. The most 
important resonance is the one at ££ m =0.33 MeV (^=10.945 MeV, J n =2 + , 3"). The data do 
not allow a unique spin assignment for this level; this uncertainty represents the main uncertainty 
in the spectroscopic factor and resonance strength. These parameters were used to estimate the 
resonance strengths for the not yet observed low energy resonances; the results are listed in Table 3. 

The reaction rate for 22 Ne(a, 7) 26 Mg is shown in figure 2 as a function of temperature. For 
temperatures below 0.3 GK the resonance at £^ m =0.33 MeV clearly dominates the reaction rate, 
possible contributions can also come from the two 1~ resonances at £^ m =0.538 MeV (^=11.153 
MeV) and ££ m =0.568 MeV (^=11.183 MeV) at the neutron threshold. Based on the recent re- 
analysis of the neutron capture data (Koehler 2002) these states are expected to have only small 
neutron partial widths. The reaction rate in the temperature range T>0.3 GK is determined 
by the contributions of states above 11.3 MeV, in particular by the strong level at -E^, m =0.711 
MeV (E x =l 1.326 MeV) which has been measured by Wolke et al. (1989). Also indicated is the 
uncertainty range for the 22 Ne(a, 7) 26 Mg rate. It is obvious from this plot that the experimen- 
tal uncertainties associated with the low energy resonances dominate the uncertainty in the low 
temperature range. 
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3.3. The comparison of the reaction channels 

The temperature dependence of the two 22 Ne+ a reaction rates differ considerably at lower 
temperatures because the neutron channel opens only at higher energies. Figure 2 demonstrates 
that at low temperature T<0.2GK the (0,7) reaction channel dominates over the competing (a,n) 
reaction process. However, there are still considerable uncertainties associated with the low tem- 
perature range of the reaction rates as demonstrated in figure 3. Shown are the upper and lower 
limits of the reaction rate normalized to the recommended value and clearly demonstrates that be- 
low T<0.3 GK the uncertainty range rapidly increases towards lower temperatures to one order of 
magnitude for the 22 Ne(a, n) 25 Mg rate and to nearly two orders of magnitude for the 22 Ne(a, 7) 26 Mg 
rate. This affects the reliability for nucleosynthesis predictions in rapidly changing low temperature 
environments. Figure 4 shows the ratios of 22 Ne(a, n) 25 Mg rate and 22 Ne(a, 7) 26 Mg rate. A com- 
parison between the upper and lower limits of the respective reaction rates shows the uncertainty 
range for this ratio. At higher temperatures the 22 Ne(a, n) 25 Mg reaction clearly dominates while 
towards lower temperatures the 22 Ne(a, 7) 26 Mg reaction will be far stronger than the competing 
22 Ne(a, n) 25 Mg channel. However, the present data are not stringent enough to determine the exact 
temperature for this cross over but can be limited to a temperature range between 0.15 and 0.3 
GK. Further measurements of low temperature resonances in both reaction channels are therefore 
necessary to improve on the present uncertainties. 

4. Results 

We now present the results of the stellar models. Unless otherwise stated, we discuss the 
models computed using the Vassiliadis and Wood (1993) mass-loss rates. 

4.1. Structural properties and the CNO nuclei 

In this section we present some of the structural properties of the stellar models used for the 
present study in detail, including a brief discussion of the evolution of the CNO nuclei during the 
TP-AGB phase. In table 4 we first include the initial mass and metallicity and the mass-loss law 
adopted for the structure calculation. VW93 refers to models calculated with Vassiliadis and Wood 
(1993) mass loss whilst R75 refers to models calculated with Reimers (1975) with the parameter 
r/R = 3.5. Furthermore, we present the core mass at the beginning of the TP-AGB, the mass of 
the envelope when HBB is shutoff and the mass at the end of the calculation, the total number of 
TPs, the number of TDU episodes, the total amount of matter dredged into the envelope Mj°* d , the 
maximum temperature at the base of the convective envelope, the maximum temperature in the He- 
shell during a TP and the maximum dredge- up efficiency, A max . In table 5 we present some surface 
abundances results including the final C+N+O value divided by the initial (see table caption for 
more details). The dredge-up efficiency is defined by A = AMd rc dg C /AMh, where AMd re dge is the 
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amount of H-exhausted core matter mixed into the envelope and AM^ is the amount by which the 
core mass grew during the previous interpulse period. For all the models considered here, A > 0.9, 
is reached after a small number (~ 6) of pulses, similar to behavior observed by Stancliffe et al. 
(2005) and Herwig (2004a) for the same mass. 

The VW93 mass-loss rate depends on the pulsation period until the period reaches 500 days, 
after which the luminosity-driven superwind phase begins. The pulsation period (Vassiliadis and 
Wood 1993) depends in turn on the radius and total mass and because the lower metallicity models 
are more compact they take longer (in time) to reach the start of the superwind phase and hence 
experience more TPs than the solar metallicity models. In comparison, the R75 massdoss rate is 
proportional to both the radius and the luminosity and as we can see from table 4 at the lowest 
metallicity we consider (Z = 10 -4 ) the large luminosity at this composition results in fewer TPs 
than using VW93 mass loss, the reverse of the behavior observed at Z = 0.02. 

For all models except the Z = 0.02 cases, the temperature at the base of the envelope, Tbce> 
increases quickly at the beginning of the TP-AGB and exceeds 50 x 10 6 K by about the 8 th pulse. 
The envelope mass after which HBB is shutoff decreases with a decrease in metallicity, see column 5 
in table 4. For the solar composition VW93 (R75) models, Tbce reaches 50 x 10 6 K after about 
the 9 th TP (12 th TP) and decreases quickly once the envelope mass drops below ss 2.6M , which 
is when HBB is shutoff. After this, dredge- up continues increasing the C/O ratio above unity in 
some cases (see table 5). 

For the remainder of this section we summarize the nucleosynthesis of the CNO nuclei, because 
these species are not affected by varying the 22 Ne +a reaction rates. The contribution of many 
TPs with efficient dredge-up means that there is a considerable increase in the amount of 12 C that 
can be converted to primary 14 N by HBB (note also the low final 12 C/ 13 C ratios shown in table 5). 
This increase can be best appreciated by comparing the final envelope C+N+O abundance (divided 
by the initial) for each model in table 5. In the lowest Z models 14 N is the most abundant isotope, 
and an increase in Z results in decreasing levels of enrichment, where the final [ 14 N/Fe] 3 values are 
3.9, 1.9, 1.5 and 0.6 in the Z = 0.0001,0.004,0.008 and 0.02 model, respectively. [C/N] ratios have 
been observed in metal-poor stars in the halo, with the finding that —0.5 < [C/N] < 1.5 at carbon 
abundances around [C/Fe] ~ 2 (Johnson et al. 2005). Our models overproduce N compared to C, 
with [C/N] < — 1 in all metal-poor 5M computations. 

The solar-metallicity VW93 model does not become carbon rich, whereas the R75 model 
becomes a carbon star after the 28 th TP. From table 4 we see that the reason for this is more TDU 
episodes, resulting in about a factor of two more matter dredged into the envelope. The 5M Q , 
Z = 0.008 model does not become a carbon star, whereas the Z = 0.004 model does. The final 
[ 12 C/Fe] values for both cases are 0.16 and 0.70, compared to 0.09 and 0.16 initially. The Z = 0.004 



assuming solar values from Anders & Grevesse (1989), as used in the Z — 0.02 models. We use the standard 
logarithmic notation [-X/Y] = log 10 (X/Y)» — log lo (X/Y)0. 
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model has a longer TP-AGB phase (0.59 Myr compared to 0.50 Myr), leading to more 16 destroyed 
by HBB, even though HBB temperatures are similar in both models. The final [ 16 0/Fe] values are 
0.02 and 0.04, compared to the initial enhanced values of 0.16 and 0.24, respectively. The 5M Q , 
Z = 0.0001 models (with VW93 and R75 mass loss) become carbon stars quickly, after the second 
TP. This is a consequence of deep third dredge-up and a low initial 16 envelope abundance. The 
final [ 12 C/Fe] values are 2.4 and 2.6, in the VW93 and R75 models. The envelope 16 abundance 
is destroyed during the interpulse by HBB at temperatures ~ 90 x 10 6 K whereas the many efficient 
TDUs result in an overall increase, leading to a final [ 16 0/Fe] ~ 1 in both cases. This behavior was 
also observed by Herwig (2004a) for the equivalent mass and metallicity. 



4.2. Surface abundance evolution of the Mg isotopes 

The envelope abundances of the Mg isotopes are first modified by the operation of the SDU, 
and the magnitude of such modifications increases with a decrease in the metallicity or an increase 
increase in the mass. For the 5M Q , Z = 0.02, the envelope 25 Mg abundance decreases by 4.2% 
while 26 Mg increases by 3.7% after the SDU, whereas in the 5M Q , Z = 0.004 model we observe a 
decrease of 9.4% in the 25 Mg abundance and an 8.4% increase in 26 Mg. For all cases but the 5M Q , 
Z = 0.0001 model, the 24 Mg abundance does not change. For this model, the 24 Mg abundance 
increases by 11%, whereas the 25 Mg and 26 Mg abundances decrease by 13% and 6.8% respectively. 
The results quoted above are for the calculations that use the NACRE NeNa and MgAl chain 
reaction rates. 

In figure 5 we show the evolution during the TP-AGB phase of the 25 Mg (solid line) and 26 Mg 
(dashed line) abundances at the surface of the 5M , Z = 0.02 model for four different choices of 
the 22 Ne +a reaction rates. The abundances in figure 5 are the mole fraction, Y (where the mass 
fraction X = Y A and A is the atomic mass), scaled to the total magnesium abundance, Mg = 
24 Mg+ 25 Mg+ 26 Mg. In figure 6 we show the surface abundance evolution of the Mg isotopes for 
three different initial compositions, that use the new recommended rates for the 22 Ne +a reactions. 
In table 5 we also present the Mg isotope ratios (in terms of 2 Mg: 25 Mg: 26 Mg) at the beginning 
and at the end of the TP-AGB phase for each model. Figure 6 and table 5 shows the influence of 
metallicity on the evolution of the Mg isotopes, where the effects of HBB become more significant 
at lower Z. 

From figure 5 we see the important result that using the new recommended rates for the 
22 Ne + a reactions does not result in large changes to the production of 25,26 Mg compared to the 
calculations performed using older estimates. For the solar metallicity model, the final surface Mg 
isotopic ratios are 70.8:13.5:15.6 when using NACRE compared to 72.2:12.4:15.4 when using the 
new recommended rates. The only observable change is less 25 Mg but about the same amount (or 
just slightly less) 26 Mg. These results are consistently observed in the metal-poor AGB models, 
although the trends are stronger at lower metallicity. 
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To explain the above trends, we need to examine the behavior of the new rates compared to the 
NACRE and K94 rates. We produce less 25 Mg because the older estimates for the 22 Ne(a, n) 25 Mg 
reaction are up to 40% faster at typical He-shell burning temperatures i.e. 250 x 10 6 K to 400 x 10 6 K. 
Why the production of 26 Mg is reduced or remains about the same is not immediately obvious given 
that the new recommended rate is faster than both the NACRE (by up to 26%) and K94 (by 70% 
to 100%) rates. The reason is again related to the (a, n) reaction which releases neutrons during a 
TP, some of which are captured to produce 26 Mg, a process that is most efficient in the lowest Z 
models. In the solar metallicity models the production of 26 Mg is less dependent on the neutron 
flux and we produce about the same when using the new recommended rate for the 22 Ne(a, 7) 26 Mg 
reaction. 



4.3. Stellar yields 

In the first row of table 6 we present the stellar yields of 25 Mg and 26 Mg (in M Q ) for the 
5M models of different Z, computed using the new recommended rates. In other rows we present 
the percentage difference between yields calculated with a different set of rates (as noted in the 
table) and our reference yields given in the first row 4 . In figure 7 we show the percentage difference 
between the stellar yields calculated using the new recommended rates (the reference), the new 
upper and lower limits and the NACRE recommended rates, for two different metallicities. The 
x-axis is atomic mass, and we present percentage differences for 22 Ne through to 31 P, with 26 Al 
included at 26.5. According to the definition of the percentage difference given above, a positive 
difference means that that model produced (or destroyed) more of a given species compared to our 
reference case. To compare to the magnitude of the uncertainties associated the NACRE rates, 
in figure 8 we show the percentage difference for models that use the NACRE recommended (the 
reference) and the NACRE upper and lower limits. We also vary the 25 Mg(n, 7) 26 Mg reaction rate 
by a factor of two each way (see figure 8). 

These figures show that using the new rates for the 22 Ne + a reactions results in considerably 
smaller uncertainties in the production of all isotopes between 22 Ne to 31 P compared to models 
that use the NACRE rates. For the 25 > 26 Mg isotopes, differences of the order ~400% in the case 
of the 5M , Z = 0.02 model with the NACRE upper limits are reduced to less than 30%, when 
using the new upper limits. For models of lower metallicity, the uncertainties are even smaller (at 
most ~ 20%). Varying the 25 Mg(n, 7) 26 Mg reaction rate given by Weigmann et al. (1976) by a 
factor of two each way results in differences of at most 33% in the production of 31 P, about 26% 
in the production of 25 Mg and about 15% in the production of 26 Mg for all models, regardless of 
metallicity. Note that the abundance of heavier isotopes are affected by neutron captures hence 
the dependence on these reaction rates. 



4 according to [yield(i) — yield(ref)]/yield(ref) multiplied by 100, where i is one of the combinations of reaction 
rates described in §2 
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The percentage differences reported in table 6 and shown in figures 7 and 8 illustrate that the 
computed errors decrease with a decrease in the metallicity of the model. The reason could be that 
the yields significantly increase with an decrease in Z, and hence changes to a larger number result 
in smaller differences than changes to a smaller number (in the case of the Z = 0.02 models). Also, 
the low Z models have more TPs and more efficient HBB means that the proton-capture MgAl 
chain reactions have a greater overall effect on the stellar yields. 

4.4. The effect of partial mixing zones and varying the mass-loss rate 

The inclusion of a partial mixing zone (PMZ) at the deepest penetration of the TDU will mix 
protons from the envelope into the He-intershell, producing a 13 C pocket. In the PMZ neutrons 
are liberated during the interpulse period by the neutron source reaction 13 C(a,n) 16 and the Mg 
isotopic abundances can be modified because of the activation of the chain of neutron captures 
starting from the abundant 22 Ne and proceeding through 23 Na to 24 Mg, 25 Mg and 26 Mg. Given the 
high value of the neutron capture cross section of 25 Mg, with respect to those of 24 Mg and 26 Mg 
(see Bao et al. 2000), neutron captures in the 13 C pocket produce 24 Mg and 26 Mg, while depleting 
25 Mg. The details of how the 13 C pocket forms is still unknown although various mechanisms have 
been proposed; see Lugaro et al. (2004) for a wider discussion. In nucleosynthesis studies, the 
extent of the pocket is usually set as a free parameter with typical values ~ 1 /15 th the mass of the 
He-intershell (Gallino et al. 1998; Goriely & Mowlavi 2000). In intermediate mass stars, the mass 
of the He-intershell is smaller by about an order of magnitude compared to lower mass stars and 
hence the importance of the 13 C pocket may be lessened (Gallino et al. 1998). As done in previous 
nucleosynthesis studies (Gallino et al. 1998; Goriely & Mowlavi 2000; Lugaro et al. 2004), we 
artificially include in the post-processing calculation a PMZ of constant mass M pmz = 1 x 1O _4 M0 
in the 5M , Z = 0.02 model. We chose a proton profile in which the number of protons decreases 
exponentially with the mass depth below the base of the convective envelope exactly in the same 
way as described in Lugaro et al. (2004). The mass of the He-intershell decreases with evolution 
and for the solar metallicity model, the final He-shell mass is about 0.001, so that our PMZ region 
is at maximum l/10 th of the He intershell. 

The inclusion of the PMZ produces small percentage differences (less than ~10%) for most 
species, with the exception of 31 P (28%) and the particle g (104%). Hence there are no changes 
to the 25 ' 26 Mg yields with the introduction of a PMZ of M pmz = 1 x 10~ 4 M Q . This conclusion is 
consistent with the results obtained by Lugaro et al. (1999) for the elements Si and Ti and can be 
extended to all intermediate-mass elements. There are two reasons for this. First, since elements 
lighter than iron have very small neutron capture cross sections, as much as 3 orders of magnitude 
smaller as compared to those of nuclei heavier than iron, their abundances are not strongly affected 
by the neutron flux in the 13 C pocket. In fact, the changes in the abundances of intermediate mass 
elements in the 13 C pocket are at maximum of one order of magnitude, while the abundances of 
heavy elements are produced by up to three orders of magnitudes of their initial values (see e.g. 
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Table 2 of Lugaro et al. (1999)). Second, the 13 C pocket is engulfed and diluted (of factors between 
l/20 th to l/10 th ) by the next growing convective pulse, and is thus mixed with material already 
a processed during the previous pulses, and with the ashes of the H-burning shell. During this 
dilution process the signature on intermediate mass elements of neutron captures in the 13 C pocket 
is completely lost. As for proton captures occurring in the PMZ region just after the mixing of 
protons has occurred, they are expected to destroy 25 Mg only in the upper layer of the PMZ where 
the number of protons is higher than 0.5 (see Figure 1 of Goriely & Mowlavi (2000)). Thus they 
do not produce a strong effect when this region is diluted in the He intershell, unless the proton 
profile is very different from the one assumed here. 

We also experiment with varying the mass-loss rate. The results for the Z = 0.02 case with 
R75 mass loss are shown in figure 9. In this figure we observe large differences of a few hundred 
percent for most species. For the 5M , Z = 0.0001 R75 model we use two different choices of 
the 22 Ne + a reaction rates, as described in figure 10. The R75 model has fewer TPs than the 
VW93 case which results in smaller yields, hence the negative values. We see from figure 10 that 
changing the reaction rates has only a small effect on the yields, compared to the change induced 
by varying the mass-loss rate. This was also well demonstrated by Ventura & D'Antona (2005b) for 
intermediate mass AGB stars of Z = 0.001. The species most affected by the change in mass-loss 
law are 23 Na, 24 Mg and 26 Al, because the abundances of these isotopes depend on the duration of 
the HBB phase. 



4.5. Other model uncertainties 

With regards to the reaction rates, uncertainties remain related to the proton capture rates for 
the NeNa and MgAl chains that modify the abundances of these nuclei because of HBB. Table 6 
shows that updating these proton capture rates from our standard set to those recommended by 
NACRE did not produce any major modifications on the Mg isotopes, since variations are only a 
few percent. The variations are small because the reaction rates in the MgAl chain which affect the 
Mg isotopic ratios the most are similar between NACRE and our standard choices, with differences 
typically less than 40%. A detailed study of the impact of the uncertainties of proton capture 
reactions on the results produced by hot bottom burning on the Ne, Na, Mg and Al isotopes is 
underway (Izzard, Lugaro, Karakas & Iliadis, 2005, in preparation). 

We discussed in the introduction that major uncertainties are still present in the computation 
of the stellar models, in particular in relation to the treatment of convection and of convective 
boundaries. When our results for the 5M , Z = 0.004 model, which are calculated using MLT 
with a = 1.75, are compared with those produced by a 5M , Z = 0.001 model computed using the 
"Full Spectrum of Turbulence" (FST) prescription for convective regions by Ventura & D'Antona 
(2005a), we find that the final abundances of the Mg isotopes are lower in the FST case with respect 
to our case by about 2 dex for 24 Mg, 1 dex for 25 Mg and 0.4 dex for 26 Mg (including the contribution 
of radioactive 26 Al). In particular, in our case, 25 Mg and 26 Mg are produced in similar abundances 
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by the effect of the Ne + a reactions and deep TDU, while, in the FST case, the production of 
nuclei of mass A=26 ( 26 A1 and 26 Mg) is favored with respect to that of 25 Mg by the combined 
effect of strong HBB and weak TDU (Ventura P. 2005, private communication). These different 
behaviors should be tested against data from pre-solar grains and, using GCE models, against the 
observations of the Mg isotopes in stars, and also the galactic abundance of the radioactive nucleus 
26 Al, as obtained by the galactic gamma line at 1.809 MeV (Diehl et al. 2004). 

Herwig (2004b) implements a diffusive convective overshoot scheme at all convective borders 
in his AGB models (see also Herwig et al. 1997), which results in a considerably different stellar 
structure to our models of the same mass and metallicity. The most significant differences are hotter 
TPs and deeper TDU (where A > 1), as well as the mixing of some CO core material into the flash- 
driven convective pocket. If we compare the surface 25,26 Mg mass fractions after the 14 th TP, we 
find that these isotopes are enriched by up to 5 times more in Herwig's 5M©, Z = 0.0001 model 
compared to our equivalent case with the new recommended rates. These differences are much larger 
than those introduced by the 22 Ne +a capture rates. Herwig also reports A( 26 Mg)/A( 25 Mg) > 1 
after 14 TPs, whereas after the same number of pulses we find this ratio to be less than unity. This 
is because hotter TPs result in a stronger activation of the 22 Ne(a, n) 25 Mg reaction hence more 
neutrons are released that are captured by 25 Mg to form 26 Mg. 



5. Discussion & Conclusions 

The main important result of the present study is that the reduction of the uncertainties on 
the 22 Ne + a reaction rates has allowed us to considerably reduce the uncertainties coming from 
these rates on the production of isotopes from Mg to P in AGB stars of intermediate mass. The 
uncertainties on the Mg yields are now at a level of ~ 30%, much lower than those obtained when 
using the NACRE upper or lower limits. The yields of 25 Mg and 26 Mg are 20% to 45% and 9% to 
16%, respectively, smaller with the new rates, as compared to NACRE. The uncertainties have also 
been reduced for species heavier than Mg, where for example, the uncertainty on the production 
of P in solar metallicity models is now at a level of 35%, much lower than the 400% obtained when 
using the NACRE upper limit. These results are clearly illustrated by figures 7 and 8 and in table 6. 

From our analysis it would appear that among the uncertainties related to the stellar models, 
those coming from the treatment of convection and of the mass loss are the largest. We described 
earlier that these are enormous when compared to the current uncertainties coming from the 22 Ne 
+ a reaction rates. Much effort is needed to improve the situation for AGB models, in particular 
with respect to convection, by trying to evaluate and reduce the uncertainties, perhaps by exploiting 
all the available observational constraints. 

Our new evaluation of the 22 Ne +a reaction rates will encourage much future work, as these 
rates are important for many nucleosynthetic processes and sites. The 22 Ne(a, n) 25 Mg reaction is 
an important source of neutrons both during the final evolutionary stages of massive stars, and 
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during the AGB phase of low to intermediate mass stars. It is responsible for the production of 
heavy s-process elements in these environments (Gallino et al. 1998; Rauscher et al. 2002). The 
new rate and its uncertainties have to be tested in relation to this process. Moreover, the smaller 
uncertainties of the rates presented here with respect to those in NACRE, appear to rule out the 
possibility that the production of the relatively abundant p-only isotopes of Mo and Ru could be 
related to a high value of the 22 Ne(a, n) 25 Mg rates (Costa et al. 2000). 

It is important to study the relative production of 25 Mg and 26 Mg, as we have done in table 6, 
because both spectroscopic observations and the analysis of pre-solar grains are able to separate 
these two isotopes. Also the contribution of radioactive 26 Al to the abundance of 26 Mg has to be 
carefully evaluated. In particular, one pre-solar spinel grain (Zinner et al. 2005, OC2) appears to 
bear the signature of nucleosynthesis in intermediate mass AGB stars, with excesses in both 25 Mg 
and 26 Mg. A future application of our present work will be to compare our detailed results to the 
composition of this grain, extending the study to the oxygen isotopic ratios in AGB stars and their 
uncertainties, as these are also measured in the grain. 

Acknowledgments 

The authors wish to thank the anonymous referee for many helpful comments that have im- 
proved the clarity of the paper. AIK wishes to acknowledge the Canada Foundation for Innovation 
(CFI) and the Nova Scotia Research and Innovation Trust fund (NSRIT) for partly funding com- 
putational resources used for this study. Financial support from R.G. Deupree's Canada Research 
Chair (CRC) fund is gratefully acknowledged. AIK warmly thanks Maria Lugaro, Brad Gibson and 
John Lattanzio for their hospitality, and acknowledges the Institute of Astronomy (University of 
Cambridge), the Supercomputing and Astrophysics department (Swinburne Institute of Technol- 
ogy) and the CSPA (Monash University) for travel support. We thank Paolo Ventura for sharing 
unpublished results on the Mg isotopes. 

A. Details of the reaction rates used in the reference case 

Most of the 506 reaction rates come from the REACLIB Data Tables (1991 version). The 
updates made to the proton, a and neutron capture rates are detailed in Lugaro et al. (2004) with 
one exception. We have since updated the 14 N(a, 7) 18 F reaction rate to that given by Gorres et al. 
(2000). This reaction rate set is considered our standard set. Note that in our standard reaction 
rate set, we use the 22 Ne +a rates given by Kiippeler et al. (1994) and Drotleff et al. (1993). 

We have also used NACRE rates for many of the reactions involved in the NeNa and MgAl 
chains. In Table 7 we list the reactions that we changed to those given by NACRE; we also include 
the reference for that rate used in the standard reaction rate set for comparison. CF88 refers to 
Caughlan & Fowler (1988). 
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Table 1: Level parameters for a unbound states in 26 Mg. 
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3 b 




6.5^-9. 10- 15 


3 


11.112 


0.497 


2+ 


1.73 


2577 


7.3l5 4 3 5 -10- n 


2580 


11.153 


0.538 


1- 


4.79 


14.4 


4.ll^-10- 9 


19.2 


11.163 


0.548 


2+ 


4.56 


4640 


8.4l^. 10 -io 


4650 


11.171 


0.556 


2+ 




1.44 


i.2i^-io- 9 


20 


11.183 


0.568 


l - 


3 b 


0.54 


7.61^.10-9 


3.54 


11.194 


0.579 


2+ 


3 b 





3.ili^-io- 9 


3 


11.274 


0.659 


2+ 


3.24 


540 


6.11U-10- 8 


543 


11.286 


0.671 


1- 


0.79 


1256 




1257 


11.310 


0.695 


1- 


3 b 


1.12 


5.8li;|-10- 6 


4.1 


11.326 


0.711 


l - 


3 b 


0.6 


9.4i^-10" 6 


3.6 


11.328 


0.713 


l - 


3.6 


424 


3.9±iJ-10- 5 


428 



"reference Koehler (2002) 

'estimated average value (Koehler 2002) 

Reference Giesen et al. (1994) 
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Table 2: Recommended resonance strengths as well as their upper and lower limits for the 



Ne(a,7) 26 Mg reaction derived from the parameters listed in table 1 



x (MeV) 


E^ m (MeV) 


J n 


w 7(a, 7 ) (eV) 




(eV) 


10.693 


0.078 


4+ 


1.6-10" 45 


1.2-10~ 46 


3.6-10" 44 


10.945 


0.330 


2+,3 - 


2.8-10" 15 


1.8-10~ 15 


1.0-10~ 13 


11.112 


0.497 


2+ 


2.4-10' 13 


3.4-10~ 14 


7.3-10~ 13 


11.153 


0.538 


1- 


3.0-10~ 09 


2.1-10" 10 


8.7-10~ 09 


11.163 


0.548 


2+ 


4.M0~ 12 


2.4-10~ 12 


6.7-10- 12 


11.171 


0.556 


2+ 


8.9-10~ 10 


4.5-10- 10 


1.3-10- 09 


11.183 


0.568 


1- 


1.9-10- 08 


1.3-10- 08 


2.4-10- 08 


11.194 


0.579 


2+ 


1.6-10- 08 


7.5-10- 09 


2.4-10- 08 


11.274 


0.659 


2+ 


1.8-10- 09 


1.2-10- 09 


2.5-10- 09 


11.286 


0.671 


1- 


6.4-10- 10 


4.2-10- 10 


8.7-10- 10 


11.310 


0.695 


l~ 


3.6-10~ 05 


3.2-10~ 05 


4.0-10" 05 


11.326 


0.711 


1~ 


2.3-10~ 05 


1.6-10- 05 


3.2-10~ 05 


11.328 


0.713 


l - 


9.9-10~ 07 


6.8-10~ 07 


1.3-10" 06 
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Table 3: Recommended resonance strengths as well as their upper and lower limits for the 
22 Ne(a, n) 25 Mg reaction derived from the parameters listed in Table 1. 



x (MeV) 


E^ m (MeV) 


r 


W7(a,n) ( eV ) 






11.112 


0.497 


2+ 


3.6-10- 10 


5.0-10- 11 


l.i-io- 09 


11.153 


0.538 


l - 


9.2-10~ 09 


6.3-10~ 10 


2.6-10~ 08 


11.163 


0.548 


2+ 


4.2-10" 09 


2.4-10" 09 


6.8-10~ 09 


11.171 


0.556 


2+ 


4.3-10- 10 


2.2-10~ 10 


6.5-10~ 10 


11.183 


0.568 


\~ 


3.5-10" 09 


2.3-10~ 09 


4.6-10- 09 


11.274 


0.659 


2+ 


3.0-10- 07 


l.O-lO" 07 


4.1-10- 07 


11.286 


0.671 


1- 


1.0-10~ 06 


6.6-10- 07 


1.4-10- 06 


11.310 


0.695 


1- 


4.7.IO-O6 


3.5-10~ 06 


6.0-10- 06 


11.326 


0.711 


\~ 


4.7-10" 06 


3.1-10" 06 


6.5-10~ 06 


11.328 


0.713 


\~ 


1.18-10" 04 


1.03-10" 04 


1.33-10" 04 
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Table 4: Structural properties of the 5M AGB models, see the text for details. 



z 


M 


M c (l) 
(M ) 


M HBB 

(M ) 


M f 
(M ) 


No. 
TP 


No. 
TDU 


^d?ed 

(M ) 


T-imax 
J BCE 

(10 6 K) 


T-imax 
^He 

(10 6 K) 


Amax 


0.02 


VW93 


0.861 


2.564 


1.499 


24 


22 


5.027(-2) 


64 


352 


0.961 


0.02 


R75 


0.861 


2.560 


1.802 


38 


35 


1.047(-1) 


57 


368 


0.977 


0.008 


VW93 


0.870 


1.857 


1.387 


59 


56 


1.745(-1) 


81 


366 


0.952 


0.004 


VW93 


0.888 


1.560 


0.944 


83 


80 


2.250(-l) 


85 


379 


0.959 


io~ 4 


VW93 


0.910 


1.056 


0.572 


137 


134 


3.133(-1) 


92 


380 


0.980 


io~ 4 


R75 


0.909 


1.100 


0.335 


70 


67 


1.569(-1) 


90 


383 


0.950 
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Table 5: Some surface abundance results for the 5M Q models including the final C+N+O and 
12 C/ 13 C ratio, and C/O ratio and the Mg isotope ratio ( 24 Mg: 25 Mg: 26 Mg) at the beginning and 
end of the TP-AGB phase (using the new recommended rates). 



z 


M 


CNO// 


12 C/ 13 C f 


C/O ag b 


c/o f 


Mg Ratio ag b 


Mg RatiOf 






CNO 












0.02 


VW93 


1.341 


7.828 


0.291 


0.766 


79.0:09.6:11.4 


72.2:12.4:15.4 


0.02 


R75 


1.814 


15.57 


0.291 


1.705 


79.0:09.6:11.4 


56.1:18.0:25.9 


0.008 


VW93 


2.886 


7.383 


0.239 


0.660 


82.8:07.6:09.6 


46.1:19.6:34.2 


0.004 


VW93 


5.106 


10.77 


0.204 


2.102 


85.3:06.4:8.33 


31.7:21.6:46.7 


10~ 4 


VW93 


431.1 


9.223 


0.209 


9.400 


82.2:08.2:09.6 


11.7:17.5:70.7 


10" 4 


R75 


339.4 


12.15 


0.205 


18.60 


82.2:08.2:09.6 


04.5:16.7:78.7 
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Table 6. In the first row we show the yields of 25 Mg (in roman font) and 26 Mg (in italics) in solar 
masses for different stellar models (M(M ), Z) that use the new recommended rates for the 22 Ne 
+q reactions. In other rows, we show the percentage difference between models computed using 
other estimates of the 22 Ne +a reactions and the yields from row 1 (see definition in the text). 



22 Ne(a, 7 ) 26 Mg 


22 Ne(a,n) 25 Mg 


(5,0.0001) 


(5,0.004) 


(5,0.008) 


(5,0.02) 


recommended 


recommended 


5.80(-4) 


4.84(-4) 


3.49(-4) 


6.49(-5) 






2.28(-3) 


l.ll(-3) 


6.84(-A) 


1.17(-A) 


standard™ 


standard™ 


-3.42 


29.6 


36.6 


45.6 






9.05 


7.50 


3.51 


-5.95 


K94 b 


K94 6 


22.5 


31.5 


38.2 


52.6 






9.80 


8.64 


5.68 


-2.99 


NACRE C 


NACRE C 


20.4 


27.9 


34.6 


44.8 






15.3 


16.0 


15.9 


8.82 


upper limit 


upper limit 


3.34 


8.94 


11.5 


24.7 






13.8 


17.0 


18.9 


26.7 


lower limit 


lower limit 


-8.93 


-7.95 


-9.98 


-14.7 






-14.2 


-16.2 


-17.5 


-16.6 


lower limit 


upper limit 


5.10 


11.5 


13.5 


20.7 






0.46 


0.21 


-1.46 


-4.13 


upper limit 


lower limit 


-9.95 


-8.94 


-11.2 


-16.5 






2.47 


3.82 


4.96 


7.82 



"Standard Case: Kappeler et al. (1994, K94) for the 22 Ne +a reactions; reaction rates for the NeNa and MgAl 

chains described in the Appendix 
6 K94: K94 for the 22 Nc +a reactions; NACRE recommended rates for the NeNa and MgAl chains 
C NACRE: NACRE recommended rates for the 22 Ne +a reactions and the NeNa and MgAl chains 
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Table 7. Proton capture rates changed to those given by NACRE 

reaction reference 

19 F(p, 7 ) 20 Ne CF88 

20 Ne(p, 7 ) 21 Ne CF88 

21 Ne(p, 7 ) 22 Na El Eid & Champagne (1995) 

22 Ne(p, 7 ) 23 Na a > fe El Eid & Champagne (1995) 

23 Na(p, 7) 24 Mg El Eid & Champagne (1995) 

23 Na(p, a) 20 Ne El Eid & Champagne (1995) 

24 Mg(p, 7 ) 25 A1 Powell et al. (1999) 

25 Mg(p, 7 ) 26 AF/ i Iliadis et al. (1996, 1990) 

26 Mg(p, 7 ) 27 A1 Iliadis et al. (1990) 

26 Al»(p, 7 ) 27 Si Champagne et al. (1983); Vogelaar et al. (1996) 

26 AP(p, 7 ) 27 Si CF88 

27 Al(p, 7 ) 28 Si Iliadis et al. (1990); Timmermann et al. (1988) 

27 Al(p, a) 24 Mg Timmermann et al. (1988); Champagne et al. (1988) 

a but beware the error in the NACRE analytical fit for 0.15 > Tg > 2.0 compared to the tabulated 
rate. 

b The rate given by El Eid & Champagne (1995) contains a typographical error for the first term. 
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FIGURE CAPTIONS 

Fig. 1. — The stellar yields of (a) 25 Mg and (b) 26 Mg (in M ) as a function of mass and metallicity, 
from Karakas & Lattanzio (2003). Also included are the yields from the Z = 0.0001 AGB models. 
The reaction rates used for these models are our standard choice as outlined in the text. We define 
the stellar yield to be yt = f T [X(k) — Xo(k)](dM / dt)dt, where r is the stellar lifetime, X(k) is the 
current mass fraction, Xo(k) the initial mass fraction, and dM/dt is current mass- loss rate. 

Fig. 2. — Reaction rates for 22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg as a function of temperature. The 
solid lines represent the uncertainty range in the reaction rate as discussed in the text. 

Fig. 3. — Ratio of the upper (lower) limits of the 22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg reaction rates 
and the recommended values as a function of temperature. 

Fig. 4. — Ratio of the reaction rates of 22 Ne(a, n) 25 Mg and 22 Ne(a, 7) 26 Mg as a function of tem- 
perature. Shown are the ratio between the recommended values and the ratio between the upper 
limit of 22 Ne(a, n) 25 Mg and the lower limit of 22 Ne(a, 7) 26 Mg as upper limit for the overall uncer- 
tainty ratio as well as the ratio between the lower limit of 22 Ne(a, n) 25 Mg and the upper limit of 
22 Ne(a, 7) 26 Mg as a lower limit for the ratio. 

Fig. 5. — The evolution of the 25 Mg (solid line) and 26 Mg (dashed line) abundances at the surface 
of the 5M , Z = 0.02 VW93 model for four different choices of the 22 Ne + a reaction rates. 

Fig. 6. — The surface abundance (in log mole fraction) of the Mg ( 24 Mg dotted line, 25 ' 26 Mg as in 
Fig. 5) isotopes as a function of time during the TP- AGB for models of different metallicity, using 
the new recommended rates. 

Fig. 7. — Percentage difference between stellar yields computed with the new recommended rates 
(the reference) and the new upper and lower limits and the NACRE recommended rates for the 
22 Ne + a reactions. We show results for the 5M , Z = 0.02 (a) and for Z = 0.0001 model (b). 

Fig. 8. — Same as figure 7, but now the reference yields were computed with the NACRE 22 Ne +a 
recommended rates and we show percentage differences for models using the NACRE upper and 
lower limits. We also show results where the 25 Mg(n, 7) 26 Mg reaction rate was varied by a factor 
of two each way. 

Fig. 9. — The percentage difference between the yields from the 5M , Z = 0.02 model with R75 
mass loss. 

Fig. 10. — The percentage difference between the yields from the 5M , Z = 0.0001 model with 
VW93 mass loss and two models with R75 mass loss. 
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